Supersymmetric hybrid inflation with a hght scalar 
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Light scalar fields present during inflation can lead to interesting observable signatures, especially 
in models with non-equilibrium reheating dynamics. We study a supersymmetric hybrid inflation 
model with a third scalar whose lightness is protected by symmetry, using analytical and rmmerical 
techniques, and demonstrate that the amplitude, the spectral index and the non-Gaussianity pa- 
rameter /nl of the primordial curvature perturbation are within the Planck observational bounds 
for suitable parameter values. 



I. INTRODUCTION 

The recently announced measurements of the cos- 
mic microwave background (CMB) anisotropies by the 
Planck satellite [l|, |2| provided even more evidence for 
inflation, a period of accelerating expansion in the early 
universe. This acceleration is believed to have been 
caused by a scalar field called inflaton, but so far obser- 
vations have shed little light on the detailed properties 
of this field. 

In particular, the inflaton should ultimately be incor- 
porated as part of the same quantum field theory that 
describes particle physics at high energies beyond the 
Standard Model. One of the most promising attempts 
to achieve this is hybrid inflation [3| embedded in super- 
symmetric (SUSY) theories such as F-term versions in 
iV = 1 SUSY, D-term versions in SUGRA, P-term ver- 
sions in A^ = 2 SUSY, D-brane versions and many more 
(see 0, [3 for a few representative examples) . 

Recent work has also demonstrated that inflationary 
models with light scalar fields, besides the inflaton field, 
can have interesting and distinctive observable signa- 
tures. The simplest example is the generation of gen- 
eration of curvature perturbations through the curvaton 
mechanism [gj] . If the light scalar is coupled to the infla- 
ton and reheating at the end of inflation involves non- 
equilibrium processes the observational signatures are 
richer, including anisotropies in the gravitational wave 
background [7| and highly non-Gaussian contributions to 
the curvature perturbation [Sj. If the perturbations are 
subdominant, they would generally lead to a small non- 
Gaussianity parameter /nl on large scales |9|| and would 
therefore be compatible with the Planck data. 

The effects of non-equilibrium dynamics of light scalars 
have been studied mostly in the context of preheating 
with a parametric resonance, but analogous observable 
signatures should also be produced during tachyonic pre- 
heating |10l . |ll | in hybrid inflation models. Depending 
on the details of the theory, tachyonic preheating can in- 
volve highly non-trivial non-equilibrium phenomena such 
as formation of topological defects, Q-balls or oscfllons, 
which can all be influenced by the light scalar field and 
therefore lead to observable signatures. 

The simplest hybrid inflation model consists of two 



fields, the infiaton cj) and the waterfall field x- When in- 
flation ends rapidly, which is required for non-equilibrium 
processes, the waterfall held has to be heavy, and there- 
fore it would have no effect on cosmological scales. In 
bosonic models it is possible to add another scalar held 
by hand, but if the new held is coupled to the two other 
flelds, it is generally not light. It is therefore interesting 
to consider the scenario in the context of SUSY hybrid 
inflation in which case the lightness of the field is pro- 
tected by supersymmetry. 

The aim of this paper is to investigate a simple mod- 
ification of the standard F-term SUSY hybrid infiation 
model which has a light scalar field. We analyze an infia- 
tionary scenario in which the waterfall field is heavy and 
the symmetry breaking phase transition is fast, in accor- 
dance with standard hybrid infiation. We study the pri- 
mordial perturbations generated on super-horizon scales 
and find that the model can produce the observed am- 
plitude and spectral index. 



II. THE MODEL 

Hybrid inflation models derived from SUSY can be of 
two types, F-term and D-term hybrid inflation (for a re- 
view of inflationary model building from SUSY see |l2| ) . 
Of the two, the F-term type attracts more attention be- 
cause it naturally fits with the Higgs mechanism. In its 
simplest form, F-term SUSY hybrid inflation is given by 
the supcrpotential [J] 
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where $ is a gauge singlet containing the inflaton and 
A, A is a conjugate pair of superfields transforming as 
non-trivial representations of a gauge group. It is the 
most general form of supcrpotential consistent with the 
R-synimetry, under which W — >■ e^'^W, $ — > e^'^^ and 
AA is invariant. 

The scalar potential is given by [IJ] 
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where the sum is taken over all the superfields S*". For 
the superpotential (H]), this gives 
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It has a single global minimum, 
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which preserves SUSY. 

At tree level, there is no term in Eq. ([3|) that would 
drive $ to zero. It is no longer the case when radiative 
corrections are taken into account. For |$| > SUSY is 
broken, meaning that one-loop corrections to the poten- 
tial do not cancel. They are given by [4'] 



A^-E 



647r2 



■M/ln 



Ml 

A2 



(5) 



The scalar sector of the potential is given by 



Vs = aHxX-^-\ +m2|$|2 + 2a2|ci>|2|x|2 
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Let us expand the field $ into its real and imaginary 
components: $ = 4=(0 + icr). The potential ([7]) becomes 

+ ^x'{4>'+<r') + ^Hx'-v'), (8) 

where, without loss of generality, we constrained X = 
X — x/V^: X being a real scalar field. Note that the 
fields (j) and a share the same mass, whereas % does not. 
The potential (O has three degenerate global minima 
with y = 0: a symmetry-preserving vacuum (SPV) at 



where F = 0,1 respectively for a boson or a fermion and 
Mi is the effective mass measured at a scale A. These 
corrections lift the potential in the $ direction, driving 
the system towards the global minimun. 

This model is appealing for a series of reasons, perhaps 
the main one being the possibility to embed it in the 
particle physics framework. Furthermore, the potential 
([3]) has only two parameters, a and v. 

The dynamics of the phase transition depend on the 
hierarchy between the masses of the two fields $ and X, 
with behaviour interpolating from a fast phase transi- 
tion that happens in less than one e-fold to transitions 
such that considerable expansion happens after the crit- 
ical point, predicting a red spectral index [13|. An alter- 
native regime consists of an infiationary phase dominated 
by the mass term of $ [iJl. This scenario is unlikely to 
be realised in the context of SUSY, since $'s lightness 
is protected by supersymmetry. However, considered in 
this regime, hybrid infiation generates a red spectral in- 
dex and can lead to a viable cosmology. 

Let us now modify the superpotential to introduce a 
third scalar field which is light at horizon crossing. The 
superpotential is 
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where $ is a gauge singlet. It is convenient to parame- 
terise the extra term using the dimensionless combination 
c = m/av. 

The Kahler potential is supposed to give the canonical 
kinetic terms plus negligible terms. Note that in Eq. © 
$ breaks the R-symmetry. This is not a problem because 
no observation tells us that this is a fundamental sym- 
metry of nature, and because the symmetry is restored 
as c —^ 0, it is natural to have c <C 1. 



V2r, 



X = ^ = 0, 



and two symmetry-breaking vacua (SBV) at 
= o- = 0, X = ±v- 



(9) 



(10) 



We are interested in infiationary trajectories that end 
in the SBV, with a fast phase transition driven by the 
waterfall field x- Its mass is given by m^ — —op-y^ + 
c?{i^ + 0-2) ^ \f2ara(^ and it becomes tachyonic inside 
the circle in the — ct plane defined by 
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The possible infiationary trajectories are shown in 
Figs. [1] and m The key fact is that, at large field values, 
the global attractor of the system is the SPV. Only tra- 
jectories that, on their way to the SPV, hit the tachyonic 
region m^ < 0, can potentially reach the SBV, as can be 
seen in Fig. [TJ However, hitting the instability region is 
not sufficient to guarantee that the system ends up in the 
SBV: the waterfall field x still has to roll sufficiently far 
down the ridge while in the tachyonic region. Otherwise 
the trajectory simply crosses the instability region and 
ends up in the SPV. This means the waterfall field has 
to be heavy: av ^ m or, equivalently c ^ 1. For such 
parameters, the waterfall is fast, happening in much less 
than one e-fold. 

It is easy to see from Fig. [T] that trajectories rolling 
from positive 4> end up in the SPV. This constrains the 
interesting trajectories to those with negative initial val- 
ues, (/> < 0. Fig. [2] strengthens the constraint: Only tra- 
jectories within an angle 7 from the tr = trajectory in 
the (j)— a plane hit the instability region. Assuming that 
the probability distribution of the initial conditions along 
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FIG. 1: Contour plot of the potential ^(0, x, 0). We choose 
m = V = TUpi and a = 1 for drawing purposes, although these 
parameters are unrealistic. The arrows represent possible in- 
flationary trajectories. If inflation starts at some 0i > 0, then 
the trajectory hits the SPV. On the other hand, if ^i < 0, the 
trajectory may end in the SBV. 



a circle centred at the SPV of large radius in the (p — a 
plane is uniform, the likelihood of a trajectory hitting the 
instability region is given by 
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In order to avoid fine-tuning of initial conditions, Psbv 
should not be too small, and therefore we will focus on 
the largest possible ratio c that leads to the desired phe- 
nomenology. As we shall see, this ratio turns out to be 
c sa 10~^, which gives a likelihood Psbv ~ lO^'^, and 
therefore fine tuning cannot be completely avoided. Of 
course, as usual, it may be possible to justify the fine 
tuning by anthropic reasoning. 



III. PERTURBATIONS 

To compare the model with observations, we compute 
the spectral index, the magnitude of /nl and the am- 
plitude Ar of the power spectrum using the 6N formal- 
ism |15lll6l . ll9l [. implemented with numerical and analytic 
techniques. 

The SN formalism relies upon the separate universe 
approximation, according to which points in the uni- 
verse that are causally disconnected evolve like a 'sep- 
arate universe', with dynamics following from the local 
energy density and pressure, and obeying the field equa- 
tions of a Friedmann-Robertson- Walker (FRW) universe 
(see for example [2Q|j121|). The SN formalism then notes 
that the uniform density curvature perturbation between 
two such points, Cj is given by the difference in e-folds 
of expansion between them, from some initial shared flat 
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FIG. 2: Contour plot of the potential V{(j),0,a). We choose 
m — V = TTipi and a — 1 for drawing purposes, although these 
parameters are unrealistic. The arrows represent possible in- 
flationary trajectories and the circle corresponds to m^ = 0. 
Trajectories at an angle greater than 7 to the a = line don't 
hit the m,^ < region and the SBV. 



hypersurface labelled *, to some final shared hypersurface 
of uniform density labelled /, 



C = SNi. 



(13) 



We use * to label the flat hypersurface, since in order to 
make observational predictions for an inflationary model, 
this surface must be deflned at the time when observable 
scales exited the cosmological horizon. Moreover, the 
final hypersurface / must be taken at some much later 
time when the dynamics are adiabatic and (^ is conserved. 
This can happen, for example, long after inflation ends, 
when the dynamics are dominated by a single fluid. 

When slow-roll is a good approximation at horizon 
exit, C is completely determined by the field perturba- 
tions Sip" at that time. 



C^SNiiSip^). 



(14) 



The field perturbations are extremely close to Gaus- 
sian [l3|, and if their amplitude is sufficiently small, 
the statistics of the curvature perturbation can be de- 
termined in a simple manner by Taylor expanding 
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where here and from here on we employ the summation 
convention, and subsequently we will employ the notation 
in common use, Na = dN; /dip'^. 



One then finds that the amphtude of the power spec- 
trum is given by 



and the spectral index Ug by [Ig, [l^ 

The non-Gaussianity of the perturbation is characterised 
by the reduced bispectrum /nl, which is given by ilS] 

5 No^NpNap 



(16) 



(17) 
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The current observational bounds on the observables 
at 68% confidence level are [ij 



10^ X Al = 2.19 ±0.06. 



(19) 



for the amplitude of the perturbations, 

71, = 0.9603 ±0.0073, (20) 

for the spectral index and J2| 

/nl - 2.7 ±5.8, (21) 

for the reduced bispectrum. 

IV. ANALYTICAL ESTIMATES 

In Section |V] we will calculate the observables Aq, Ug 
and /nl numerically, but in order to develop some physi- 
cal intuition about the dynamics, it is instructive to con- 
sider a simple analytical approximation first. 

The total expansion from the initial hypersurface * to 
the final uniform-density hypersurface / consists of three 
separate stages. Initially x = ^-^d both (j) < and a 
slow-roll towards the minimum of their quadratic poten- 
tial, until the trajectory hits the critical surface region 
where x becomes tachyonic, i.e., m?, becomes negative. 
We denote the amount of expansion during this first stage 
by iVcrit- It is followed by a short period of tachyonic 
growth of X, which lasts until backreaction makes the 
dynamics non-linear, and during which the universe ex- 
pands by Attach e-foldings. We assume that after this, 
the fields equilibrate instantaneously, and the universe 
expands with the equation of state of radiation to the 
final hypersurface. This expansion is denoted by A'rad- 
The total amount of expansion is therefore 



iV((/)*,CT,) = 

A'^crit (0* ,cr*) + Ntach ((/)*, Cr* ) + iVrad ( 



i= •) '^* 



(22) 



Let us first compute the number of e- folds A'crit(o'*) 
from a flat hypersurface at horizon crossing to the critical 
surface defined by rn^ — 0, which corresponds to 



= 0crit(o') 
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During this period, the fields are slowly rolling, and the 
equations of motion are therefore 
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Solving for N at the instability point, we find 

^crit(0*,</'crit,O-crit) 



(24) 
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+a2(0, - 3</>crit)mw'^ - V2a^v^^ . 

(26) 

This is the amount of infiation happening before the tran- 
sition starts and it depends on the value of cr, through 

CTcrit(f*) and (/'crit(o'crit)- 

Next, we consider the expansion Attach during the lin- 
ear tachyonic growth of x- Assuming that and x ai'e 
homogeneous, we can approximate 



0(t) W 0crit - (pt. 
a{t) « CTcrit — 0-t. 



(27) 



The time dependence of the fiuctuations x(i)k), up to a 
short time after the transition, is well approximated by 
the Minkowski linearized equation of motion 



2 

I (2a^(/)ciit + viarnjcj) + 2Q;^crcritO') t — k^ 



X(i,k) 

X(i,k). 
(28) 



Each mode k starts growing when ml. = cu^t — k'^ becomes 
negative, that is at tk = k'^/u!^. 

In a quantum field theory, Eq. (j28p is valid as an oper- 
ator equation. The initial state is the vacuum and at tree 
level it is completely described by the two-point functions 
of the fields, 

(X*(k)x(k')) = ^^i2nfS'{k-k'), 
(^*(k)7r(k')) = M(2^)3^3(k_k'), (29) 



where tt — StX- 



Following Ref. [ll| we find that, at late times, the fluc- 
tuations of X go as 
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We are interested in the time ibr at which the lin- 
ear approximation (j28l) fails and the dynamics become 
non-linear. We approximate that this happens when 

Using Eq. (|30p . the asymptotic late-time solution gives 
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(32) 



Finally, we assume that from ibr onwards, the universe 
is radiation dominated. Therefore the amount expansion 
from tbr to the final hypersurface with constant energy 
density pf is 



N, 



rad 



iln^ 
4 Pf 



where pbr = p(ibr)- From this we find 
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(33) 
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In the limit c <C v/cj)^,, we find the simple relations 



1 — 16c m-pj/w and /k 



7nl = 10c^mp[/3u^. The oppo- 
site regime, when (/)* dominates the denominator, gives 
ris = 1 - 8to2j/02 and /nl = 5m2j/(3(/)2). Therefore, 
the presence of a third scalar field which is light at hori- 
zon crossing generates a red spectral index in this model, 
in contrast to the original hybrid inflation scenario. In 
the limit in which the potential ([S]) is vacuum dominated 
at horizon crossing, this is due to the linear term in 0, 
that implies e* ~ ry* and rig < 1. As we approach 0* 
domination, the spectral index decreases. 

Planck results for the amplitude of the primordial cur- 
vature perturbation Eq. (|19p constrain the model's pa- 
rameter space. During the slow-roll phase, the amplitude 
is the only observable that depends on rescalings of the 
potential, meaning we can fix it to the observed value 
a posteriori. However, once inside the critical surface 
Eq. ITTI) . the velocity of the transition depends on the 
scale of the potential, implying that observables such as 
/nl and rig are not invariant under potential rescalings. 
Fortunately, their dependence on tJcrit does not change 
qualitatively as we vary a, and they tend to decrease 
as trajectories approach the edge of the critical surface 
Eq. (fTT|) . as can be seen in Figs. [6] and [71 The drop can 
be understood via the following physical picture. Tra- 
jectories that hit the instability close to the edge of the 
critical surface gather more kinetic energy than the other 
trajectories, slowing the expansion rate and generating a 
negative SN. This explains the decrease in /nl and the 
spectral index, as a negative SN creates a negative con- 
tribution to both observables. 



We are now ready to compute the observables of 
Eqs. (I16ll7ll8p . First, we need to calculate the deriva- 
tives Na. From Eq. (j24p we know i7*(CTcrit) and from 
Eq. ([22]) we know A^(crcrit)- Thus 
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(35) 



The analytical expressions for the observables are long 
and complicated. However, there exists a limit in which 
they become simple enough to be presented here. In- 
deed, for trajectories such that cTcrit = 0, both the first 
and second derivative of Ni are dominated by ip and the 
observables are given by 
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16c^r7ip[ 
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and 
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V. NUMERICAL RESULTS 

We have to keep in mind that the analytical results 
rely on the linear approximation (l?7|) and are therefore at 
best indicative. Nevertheless, we can use them to identify 
suitable parameters for detailed numerical study. 

In order to calculate the perturbations numerically, we 
solved the homogeneous equation of motion for the three 
fields. 



^o 



3Hipa 



dV 

dipa 



0, 



(38) 



where (pa = (</>, cr, x)- Because we are ignoring inhomoge- 
neous fluctuations, our calculations do not capture non- 
equilibrium phenomena that take place at the end of in- 
flation, which could lead to interesting effects. 

In order to solve the equations ((38)) . we have to choose 
initial conditions for the fields. Because x is heavy ini- 
tially, i.e. m^ > H, it does not have superhorizon fluc- 
tuations. Therefore inflation happens in the <f>— a plane. 
We choose the initial conditions at horizon crossing (j)* 
and (J*, and solve the equations until the trajectory os- 
cillates around the SBV. After a few oscillations, we stop 
the simulation and record the final energy density pf. 
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FIG. 3: Us as a function of c = m/{av). The continuous 
line represents the analytical result Eq. ([36} and the dots 
are numerical data points. The other parameters are fixed 
as V = 10~^mpi, (Tcrit = 10~''mpi and a is chosen to satisfy 
the Planck amplitude constraint. Although the lowest data 
point, Tis — 0.972 and c = 10"'^, is still outside the Planck 
68% confidence limit range, the spectral index decreases as 
the infiaton S becomes heavier. 



Then, we vary the initial conditions 0* and a^ and re- 
peat the simulation, ending it when the energy density 
equals pf. We repeat the process until we have enough 
points to compute the derivatives of N{(j)^:,a^,) needed to 
evaluate the observables Aq, Ug and /nl [Eqs. ([TG]) . (fT7|) 
and (|18|) respectively]. At this point, we compare the nu- 
merical Aq to the Planck constraint (obsAs) and rescale 
a as required. Finally, we repeat the process with the 
correct a, obtaining a data point. 

In practice, x ^-Iso has to be given small non-zero ini- 
tial value because x = is a classically stable trajec- 
tory protected by symmetry. Therefore if x = initially, 
the system never reaches the SBV. In reality, the sym- 
metry breaking takes place because of quantum fluctua- 
tions, but in the absence of inhomogeneous fluctuations, 
Eq. ([55]) fails to capture this. However, since our aim is 
to simply study the curvature perturbation generated on 
super-horizon scales by the light scalars (f) and ct, we only 
need to know how the time necessary for the phase tran- 
sition to the SBV depends on the field values at horizon 
crossing (j)^ and a^ . The dependence can be extracted by 
choosing a fixed initial value < Xi ^ w and varying the 
initial conditions at horizon crossing in the (j) — a plane, 
giving the function N^cjy^^a^:). 

The numerical results confirm the validity of the an- 
alytic estimates Eqs. ([55)) and (|57)) at small cTcrit, as is 
shown in Figs. |3] and S) Furtmermore, requiring the spec- 
tral index to be within the 95% confidence limit leads to 
the constraint c>5xl0~''forw = lO^^mpi. For fixed 
c, considering smaller values of v drives the spectral in- 
dex closer to unity, as can be seen from Eq. ([36]). On 
the other hand, larger values lead to smaller ng, increas- 
ing the agreement with observations. However, as can be 
seen from Eq. (j24p . increasing i; implies that (/) < rolls 
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FIG. 4: /nl as a function of c = m/{av). The continuous 
line represents the analytical result Eq. ()37p and the dots are 
numerical data points. The other parameters are fixed as 
V — 10~'^mpi, (Tcrit ~ 10~*r7ipi and a is chosen to satisfy the 
Planck amplitude constraint. 
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FIG. 5: a as a function of (Jcrit after the rescaling necessary 
to fit the observed amplitude of perturbations. The other 
parameters are fixed as v = IO~^m,pi and m = 10~^aii. The 
region to the left of the vertical dashed line leads to non- 
Gaussianity parameter /nl within the Planck 68% confidence 
limit range. 



faster through the critical surface, decreasing the cross 
section of field values ending up in the SBV. 

Eq. ([T2|) highlights that as c increases, so does the like- 
lihood of a random trajectory hitting the critical surface. 
However, increasing c has another undesired effect: the 
light scalars and a roll faster compared to x and the 
transition to the SBV might not have time to happen if 
too close to the edge of the critical surface. Indeed, we 
found that for c > 10""^, only trajectories significantly 
far from the edge hit the SBV. 

From Fig. [3| we know that for c < 10^'^ the pre- 
dicted spectral index is in conflict with observations 
([20)) in the limit that cr ^ 0. For c = 10"^, we find 
rig ~ 0.972, which is also excluded by the 68% confidence 
limit bounds, but is well within the 95% bracket. Fur- 
thermore, Fig. ID shows that /nl grows proportionally to 
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FIG. 6: Tis as a function of (Jcrit, numerical results (dashed 
line) and analytical results (continuous line). The other pa- 
rameters are fixed as a = 10~^, v — 10~^7npi and m — 
10~^aw. The region to the left of the vertical dashed line 
leads to non-Gaussianity parameter /nl within the Planck 
68% confidence limit range. 
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FIG. 7: /nl as a function of Ucrit, numerical results (dashed 
line) and analytical results (continuous line). The other pa- 
rameters are fixed as a = 10"^, v — 10~^mpi and m — 
10~^aw. The region to the left of the vertical dashed line 
leads to non-Gaussianity parameter /nl within the Planck 
68% confidence limit range. 



Considering these constraints we choose v — IQ ^rn-pi 
and c == lO^'' for more detailed study. 

Fig. [5] shows the rescaled dimensionless parameter a 
as a function of (Tcrit- It is interesting to note that it is 
at least of order 10^^, and of order 10^^ for trajectories 
leading to a viable cosmology. 

In Figs. [B] and [7] the spectral index rig and the non- 
Gaussianity parameter /nl are plotted as a function of 
fjcrit- The dots represent numerical data points whereas 
the continuous line shows the analytical results. We can 
see that the analytical approximation is only valid for 
small CTcrit- For c = 10""^, we find that the universe 
reaches the SBV at the relatively low probability Psbv ~ 
4 X 10-4. 

Requiring the spectral index to be in the observation- 



ally compatible range at 95% confidence level constrains 
the viable range of field values to |crcrit| ^ 8.3 x IQ-'^mpi. 
The non-Gaussianity parameter /nl increases the con- 
straint on the initial conditions, giving |crcrit| ^ 8 x 
lO^'^TOpi for predictions within the 68% confidence limit 
bracket. 

Note that for tJcrit > 8.3 x 10-"^?7ipi the spectral index 
drops sharply, as can be seen in Fig. [6l As discussed 
at the end of Section IIVI the sharp drop of the spectral 
index arises from the negative contribution to 5N gen- 
erated for trajectories entering the tachyonic region at 
the edge of the critical surface. These trajectories carry 
more kinetic energy inside the critical surface, leading to 
a faster transition to the SBV. 



VI. CONCLUSIONS 

In this paper we have studied a SUSY F-term version 
of hybrid inflation with three dynamically relevant scalar 
fields: the inflaton 0, the waterfall field x and a light 
scalar field cr. Softly broken R-symmetry ensures that 
the scalar a naturally has a light mass. 

The potential ([5]) has three degenerate global minima, 
one with unbroken and two with broken symmetry, as il- 
lustrated in Fig. [TJ Inflationary trajectories that end up 
in the symmetry breaking vacuum require fine tuning of 
initial conditions, which could be justified by anthropic 
arguments. For parameter values leading to observation- 
ally compatible predictions, we found that, under reason- 
able assumptions, the likelihood of an inflationary tra- 
jectory reheating in the symmetry breaking vacuum is of 
order 4 x fO~*. 

We considered the regime in which the symmetry 
breaking field x is heavy during the inflationary phase, 
i.e. m^ > H, and does not affect super-horizon observ- 
ables. We have shown that the presence of the light scalar 
a generates a red spectral ng < 1, as can be seen in 
Eq. (|36p , in agreement with observations and in contrast 
with the original hybrid inflation scenario. 

The analytical analysis carried in Section ITVl led us to 
choose parameters v = IQ-^mpi and c = m/{av) = lO"'^ 
for a detailed study. In Section |Vl we numerically com- 
puted the observables rig, /nl and A(^ for such parame- 
ters. We found that trajectories hitting the critical sur- 
face Eq. ([TT|) at values such that |CTcrit| < 8 x IQ-'^mpi 
lead to a spectral within the current 2a confidence limit 
range and to non-Gaussianity parameter /nl within the 
68% confidence limit range, as shown, respectively, in 
Figs. [HI and [T] Furthermore, for such trajectories, the 
rescaled dimensionless parameter a is of order 10"^, as 
can be seen in Fig. [5] 

It is known that the presence of light scalar fields has 
non-trivial effects on preheating dynamics 0, Q . Anal- 
ogously, light scalars should affect the non-equilibrium 
dynamics of tachyonic preheating, leading to potentially 
observational signatures. Their study requires a lattice 
simulation of the symmetry breaking phase transition. 



which is beyond the scope of this paper. 
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